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ABSTRACT 


A series of fracture tests on large-scale, pre-cracked, alumin um 

alloy panels is being carried out to examine and to characterize 

the process by which cracks propagate and link up in this 

material. Extended grips and test fixtures were specially 

designed to enable the panel specimens to be loaded in tension, 

in a 1780-kN-capacity universal testing machine. Twelve panel 

specimens, each consisting of a single sheet of bare 2024-T3 

aluminum alloy, 3988 mm high, 2286 mm wide, and 1.016 mm thick 

are being fabricated with simulated through-cracks oriented 

horizontally at mid-height. . Using existing information, a test 

matrix has been set up that explores regions of failure that are 

controlled by fracture mechanics, with additional tests near the 

boundary between plastic collapse and fracture. In addition, a 

variety of multiple site damage (MSD) configurations have been 

included to distinguish between various proposed linkage 

mechanisms. All tests but one use anti-buckling guides. At this 

writing seven specimens have been tested. Three were fabricated 

with a single central crack, three others had multiple cracks on 

each side of the central crack, and one had a single crack but no 

anti-buckling guides. Each fracture event was recorded on film, 

video, computer, magnetic tape, and occasionally optical 

microscopy. The visual showed the crack tip with a load meter in 

the field of view, using motion picture film for one tip and SVHS 

video tape for the other. The computer recorded the output of 

the testing machine load cell, the stroke, and twelve strain 

gages at 1.5 second intervals. A wideband FM magnetic tape 

recorder was used to record data from the same sources. The data 

were analyzed by two different procedures, (1) the plastic zone 

model based on the residual strength diagram, and (2) the R- 

curve. The first three tests were used to determine the basic 

material properties, and these results were then used in the 
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analysis of the two subsequent tests with MSD cracks. There is 
good agreement between measured values and results obtained from 
the model. 

KEY WORDS: aluminum, crack, fracture mechanics, multiple site 
damage, plastic zone, R-curve, toughness 


INTRODUCTION 


The aging of the commercial transport fleets around the world is 
of constant concern because of the loss of structural integrity 
through fatigue cracking. In one design approach for aircraft 
fuselages using semi-monocoque construction, circumferential 
rings or frames are intended to steer dangerous longitudinal 
cracks — if they appear — in the less threatening circumferential 
direction around the fuselage. It has been hypothesized, 
however, that in the case of aging aircraft in which multiple 
site damage, such as short fatigue cracks emanating from rivet 
holes, is present, cracks that start running longitudinally may 
continue to do so because the cracked rivet holes may provide a 
path of lesser resistance. Our research is intended to provide 
some of the information needed to better understand this crack 
propagation process. In this work we were greatly aided by the 
advice of David Broek. Also a team from NASA Langley under 
leadership of James C. Newman assisted with additional 
measurements . 

The specialized facilities and capabilities at NIST are being 
used to carry out a series of fracture tests on twelve large- 
scale, 2286-mm wide, pre-cracked, aluminum alloy panels to 
examine and to characterize the process by which cracks propagate 
and link up in this material. The tests are sponsored by the 
Federal Aviation Administration (FAA) as part of its National 
Aging Aircraft Research Program. The current tests were deemed 
necessary by the FAA because in previous work [1] only 508-mm 
wide flat panels were tested, while all other tests were 
performed on 2286-mm wide curved panels loaded by pressure, both 
with and without frames and tear straps. The results of these 
previous tests were predicted very well with a plastic zone 
model, but showed that the main mode of failure was plastic 
collapse. In the current program the wide plates failed under 
conditions closer to fracture mechanics and R-curve behavior. 

Using existing information obtained from the tests with smaller 
specimens, we set up a test matrix that explores regions of 
failure that are controlled by fracture mechanics, with 
additional tests near the boundary between plastic collapse and 
fracture. In addition, a variety of multiple site damage (MSD) 
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configurations have been included to distinguish between various 
proposed linkage mechanisms. All tests but one are planned with 
the use of anti-buckling guides. The one without anti-£“kll£g 

If buSklitgf 600 ^ 1 ’' 3 ®' 1 by HASA Langley to hel P assess the effect 


^K t ^ 1S 4- W 5 it; !'r i ? sev ? n specimens have been tested. Three were 

muittn?*» ed W1 v h 9 Slngle central crack. Three others had 
multipie cracks on each side of the central crack. One had a 

single crack with no anti-buckling guides. Each fracture event 
was recorded on film, video, computer, magnetL tape L ?he 

a ^S CC f S1 ° n ^ lly also added optical microscopy. Using flat 
sheets without stringers to stiffen the panels, these are 9 

uncomplicated tests aimed more at obtaining basic material 
properties than simulating fracture in an airplane fuselage The 

material properties sought are the basic fracture proper?!!; and 
linkage criteria for the MSD cracks. properties and 

The data were analyzed by two different procedures, ( 1 ) the 

2° ne model based on the residual strength diagram and 
(2) the R-curve. The plastic zone model is an engineering 

"5 1C V ake f the Plasticity into account by an 

fra J ture to ughness, which is less than the true 
fracture toughness of the material. The basic concept is that 
^fa_ re ® ldu ff strength in the presence of MSD depends on the 
criterion that an MSD crack will be absorbed by the main crack 

faiL bv 1 LlTa D S£ aSt T5 h° ne S meet ' S ° that the li( 3 amen t then 
£ a i~f. by c °i lapse * It has been a very good predictor of the test 
results. The R-curve uses a more fundamental approach but 

foTibT P TL r t e iIf?y |£r^ ireotly 

f £ h 5 Ge tests each contained a single crack and the 
data w ® re . used to determine the basic material 
tie f ' n ^ me1 ^ the collapse strength and the effective 
fracture toughness for the residual strength diagram, and an 
an a l y tic expression for the R-curve. These results were then 

The results a £ alysis of the two sub sequent tests with MSD cracks. 
r ?®. ® ul £f f ho ” f ood agreement between measured values and 
results obtained from the two independent analyses. 


DESCRIPTION OF TESTS 


uses 90 inch < 2286 ™> “ide panels, 
multiple site damage (MSD) , it has been named the "90 

program. The planned test matrix is shown in Table 1 


some with 
MSD" 

The 
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quantities that appear in this table are defined in Figure 2. 

The individual tests are labeled MSD— 1, HSD— 2, etc. 

Twelve panels were procured, each consisting of a single sheet of 
bare (not clad) 2024-T3 aluminum alloy, 3988 mm high, 2286 mm 
wide, and 1.016 mm thick. The specimens are being fabricated 
with simulated through-cracks oriented horizontally at mid- 
height. The simulated cracks are saw cuts, ending with the 
sharpest jeweller's saw cuts available, having a final tip radius 
of 0.076 mm. The first three tests each had a single central 
crack. Subsequent tests will also have multiple small cracks on 
each side of a larger central crack to simulate multiple site 
damage (MSD). Each MSD crack had a circular 5.6 mm diameter hole 
in its center to simulate a rivet hole, as shown schematically in 
Figure 1 . 

A review of the literature suggests that the specimens tested in 
this program are the largest structural panels that have been 
tested in tension. The great size necessitated special design 
and testing considerations in order to introduce the test loads 
uniformly along the panel widths. A whiff le-tree approach was 
ruled out by the height limitations of the 1780-kN-capacity 4- 
screw— powered universal testing machine (UTM) that was used. As 
it was, 76 mm had to be cut from the specimens reducing their 
height to 3912 mm, and only 30 to 50 mm of the testing machine's 
stroke remained at specimen failure. The machine is one of the 
largest electro-mechanical testing machines in the world. In its 
unaltered state, with power screws in all four corners of the 
1068 mm x 1524 mm testing table, loads to 448 kN can be applied 
up to 914 mm off-center, and up to full capacity at 152 mm off- 
center. 

However, to accommodate our large panels, the heads of the 
testing machine were effectively enlarged with pairs of wide 
flange structural steel beams (W8x40) , 2286 mm long, bolted 
together. The grips consisted of 2286-mm-long, thick-walled 
aluminum-alloy extrusions bolted to the steel beams. . Each end of 
the panel specimens was fastened between the grips with forty- 
five 15.875-mm high-strength steel bolts, fully tightened. The 
length of the panel between the top and bottom rows of bolts was 
3810 mm. Abrasive cloth was inserted between the specimens and 
the grips to maximize the transfer of load by shear and thereby 
avoid pin bearing failures of the thin panel material. 

The uniformity with which the load was introduced was monitored 
in the first test with 20 strain gages and in subsequent tests 
with 10 strain gages mounted on each panel about 406 mm from the 
grips at each end. These are called the far field strain gages. 
The strain distribution was measured at low loads, prior to each 
test, and, if necessary, thin metal shims were inserted between 
the steel beams and the grips in order to achieve a more uniform 
distribution. Uniformities within ±10 percent were obtained in 
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all cases. 


n2«H~? UCjCli J 9 ? uide f' consisting of four aluminum channels, were 
used to restrain out-of-plane buckling of the panel. The beams 

both «?H Ced 5°f? zontal ly about 12 mm above and below the crack on 
both sides of the specimen. In the first test a 12 mm thick felt 

smooth* 3 ?®f d between the guides and the specimen to facilitate 
smooth sliding. In subsequent tests rubber was used Test MSD-« 
^performed without the anti-buckling guides to as^ertaK ?heir 

aaaes e n? instrumented. Eight additional strain 

Msn rp^ h £ lr ? t test ) yere placed near the crack tips or 

MSD cracks. The strain gage signals were run through wide band 

instil am P lifiers - A displacement gage was 

d?f«? lled a 5 th | bottoiB . of the cross head to measure the total 
displacement. For specimens MSD-4 and MSD-5 a clipgage was also 
mounted at the middle of the central crack. The signals of all 

Y ere . ^ eco f ded by a computer. A SVHS video camera was 
mounted to view the right crack tip and a motion picture camera 
viewed the left. For tests MSD-2, 3, 4 and 6 a team from nasI 
rr 3 rv 6 f ,al f° usad high res °lution optical microscopy at the left 
hnoS^- t;LP f ^ m the backside - Test MSD-6 , which had no anti- 

? guides, used thirteen 3-element rosette strain gages. 

SSS^XToSS. 1 " the tWO 0ra0k paths and four ““ the 

The tests consisted of pulling the specimen to fracture under 

i if 5 laC ? ment i C ° n 5 ro1 * The dis Placement was generally applied at 

a ? S rT alS i ° f 2 ? t0 45 *** and held for one to four minutes at 
1 d ^ le «fi' The whole test lasted from 15 to 20 minutes. 
Link-ups to MSD cracks occurred in a fraction of a second. 

Towards the end of the test there was a large amount of crack 
growth with very small increase of load. After 50 to 100 mm of 
such tearing instability occurred and the load started to drop 
Final fracture occurred with an audible rip. P * 


DATA COLLECTION 


collection had some built-in redundancy and consisted of 
c ?®P° nen bs: manual, visual, computer, magnetic tape, 
and occasionally optical microscopy. The signals from various 
gages were recorded as a voltage and therefore prior to each test 
various calibrations were performed to be able to convert these 
voltage readings to the appropriate physical quantities. 
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Manual Recording 


During the test at each load level the load was noted . from the 
dial of the universal testing machine (UTM) . The strain at these 
loads was read with a bridge amplifier at some of the far field 
gages near the top and bottom of the specimen. For tests MSD-1 
and 7 all the far field strain gages were manually recorded and 
for MSD-6 none were manually recorded. These recordings provided 
a record of the average strain and its uniformity across the 
panel. When the final fracture occurred the fracture load was 
also recorded from the dial of the UTM. 


Visual Recording 


The visual recording consisted of showing a voltmeter with a 
crack tip and the MSD cracks (when present) in the field of view, 
using SVHS video tape for right tip and motion picture film for 
the other. The two voltmeters were connected to the UTM load 
cell and showed a voltage proportional to the load. A 
calibration was run for the voltmeters prior to the test. The 
SVHS recording could be observed live on a TV screen. Thus the 
progress of the crack growth could be monitored during the test. 
After the test the visual recordings were used to determine the 
crack extension as a function of the load. 


Computer Recording 


A personal computer was set up to collect data from various 
sources. These included the current time, the load obtained from 
the UTM load cell, the displacement, the strains from the far 
field gages that were not taken manually, the gages near the 
crack, and for tests MSD-4, 5, and 7 the clipgage. Each data set 
was taken at lh second intervals, was shown live on the monitor 
screen, and could be stored in a file. 

As mentioned the input data were read as a voltage. Therefore, 
calibrations were run prior to each test. For the load 
calibration a simulated load was generated at the UTM console 
with 22.2 kN intervals from zero up to 200 or 450 kN. The 
displacement gage was calibrated at 0.254 and 0.635 mm intervals 
over a range of 20 mm. The strain gages had a resistance of 350 
n and a gage factor of 2.135. They were calibrated using the 
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strain gage conditioning amplifiers, which also contained bridge 

°£ 3 f° n * With a switch a shunt resistance of 174?8 
* \ ° across the bridge, which corresponded to a 
simul a ted strain change of 936 microstrain. The above 
calibrations resisted in linear conversions. The clipgage worked 
on the capacitive principle and was calibrated at 0.635 mm 
intervals from zero to 15 mm. It resulted in a nonlinear 
conversion of the form a+bx+cxVx+d/x, 


Tape Recording 


A wideband FM magnetic tape recorder was used to record data from 
the same sources as the computer recording. The tape recorder 
was set to IRIG I and run at 30 inch per second, which provides a 
distortion free signal (1 dB) from DC to 20 kHz. These 
recordings are fast enough to show dynamic effects. These data 
have not yet been analyzed. 


Optical Microscopy 


MSD_2 ' 3, and 4 a team from NASA Langley collaborated 

For^it^Sn 8 ? 1 ^ 1011 2J^ al microsco PY (OM) at one crack tip. 

MSD ”f th f Y additionally recorded the signals from the 3- 
e 1 ament rosette strain gages. Their optical microscopy apparatus 

waS ° b ? erve the growing crack. A photographic technique 

Tn?^ni« P< * ♦.“S? 8 "?* tJ ?e crack tip opening angle during cirack 
initiation and stable tearing. The OM apparatus consisted of a 
computer controlled, long focal length microscope fitted to a 
video camera. The field of view of the microscope was a square 
with sides approximately 1.8 mm long. The field of view was 
centered on the crack tip and advanced as the crack grew. The 

results from this analysis will be reported elsewhere by NASA 
.Langley . J 


R— CURVE ANALYSIS 


The crack extension that occurred as the load increased was 
measured with the SVHS video tape, the movie film, and the OM 
apparatus. The results for tests MSD-l, 2, 3 and 6 are shown in 
Figure 3. These tests all had a single central crack, see Table 



1. Specimen MSD-6 had no anti-buckling guides and the data for 
this test fail below those for MSD-3, which had the same crack 
size. This shows that the lack of guides and consequent buckling 
acted as a stress intensification at the crack tips. 

We have used the R-curve concept [2] to analyze some of the data. 
First, we show here how to find the crack-extension force, G, 
from the slow stable crack growth data. This quantity is related 
to the stress-intensity factor, K, by 

G = K 2 /E (1) 

where E is Young's modulus. The stress-intensity factor is 
generally given by the generic expression 

K = jSay(Tra) (2) 

where the function f3 is used to describe the effect of the shape 
and size of the crack and the specimen. For the center crack 
specimen we then have 

G (a) = (a 2 rra/E) /8 2 (wa/W) (3) 

where a is the applied stress, a the half-crack size, and W the 
panel width. For the stress-intensity factor of a center-notch 
specimen with sharp crack tips Feddersen [3] discovered that 

/3(»ra/W) = 7sec(7ra/W) (4) 

gives an approximation that is accurate to 0.3 percent for 
2a/W<0.7. The stress is given by 

a = P/(WB) (5) 

where P is the load and B is the panel thickness. For our panels 
we have 

W = 2286 mm 
B = 1 . 016 mm 
E = 71 GPa 

Using these values and Equations (3)— (5) , the movie film and 
video tape data from Figure 3 were converted into the crack- 
extension force data shown in Figure 4. We see that for small 
crack-extensions the data from the tests with the anti -buckling 
guides collapse nicely into a master curve, but without anti- 
buckling guides the data are distinctly different. Now, 
according to the R-curve concept the crack-extension resistance, 
R, equals the crack-extension force, G, for slow stable crack 
growth with a sharp crack tip. We wish to identify those points 
for tests MSD-1, 2, and 3, shown solid in Figure 4, that qualify 
as R-curve data. The data with no crack-extension represent the 


970 



original blunt notch and so do not yet represent crack growth 
along the R-curve. For short crack-extensions up to about 75 nun 
there is slow stable crack growth, while for larger crack- 
extensions there is unstable dynamic crack growth. Initially we 
do not yet know the critical point of instability that separates 
the stable from the unstable data for each test. We therefore 
used an iterative procedure and started with an initial guess for 
the critical crack-extension. To represent the R data by a 
universal R-curve, a power law was fitted to them, giving 

R = 75.9 da 0 28 (6) 

where R is measured in kN/m and da in mm. With an analytic 
expression for the R-curve the point of instability can be found. 
The stability diagram in Figure 5 illustrates this for test MSD- 
1. A crack-extension force line or G-line is drawn and the 
control parameter is adjusted until the line is tangent to the R- 
curve when it becomes the G c -line. For simplicity we can take 
for G the expression for the infinite plate with a single crack 
(/ 3=1 ) under load or stress control 

G = a 2 7ra/E ( 7 ) 

because for cracks under 500 mm in our panels this equation does 
not differ much from the correct expression below. Equation (7) 
gives a straight line in Figure 5. Tangency is achieved by 
adjusting the stress until a critical value of o c =146 MPa is 
obtained, which then gives a critical crack-extension of da =69.3 
mm and a critical crack-extension force of G =249 kN/m. Hence, 
the data points for which 0<da<69.3 mm are stable for MSD-1 and 
those points were used for the curve fit in Figure 4. The same 
procedure was applied to the data of tests MSD-2 and 3. If it is 
found that the initial guess for da c was wrong the procedure is 
repeated until all the fitted points lie below the points of 
tangency. Figure 5 shows that for da>da the unstable points lie 
closer to the G c -line than the R-curve. 

Actually the testing machine was operated under displacement 
control so that the displacement, 6, is the control parameter. 

This quantity is related to the stress by 

a = S/( BWC t ) (8) 

where C t is the total compliance of the testing system, which can 
be decomposed into the machine, panel, and crack compliance, as 
follows 


c t = c m + H/(EBW) + 4/ (trEB) I(*a/W) ( 9 ) 

where H (150 in) is the panel height and 

I(x) = Jx/3(x)dx (10) 
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The total compliance was found for each test from the slope of 
the measured total displacement versus load curve. Using the 
slope before any crack extension occurred. Equation (9) was used 
to deduce the machine compliance. The average for the first 
three tests was C m = 1.082X10’ 5 m/kN (± 4%). The complete 
expression for the G-line under displacement control is then 
given by combining Equations (3) , (8), and (9) into 

G (a) = E$ 2 jra /9 2 (xa/W) / [ EBWC m + H + 4 W/jt I(»ra/W)] 2 (11) 


To find the critical point of instability, as was done m Figure 
5, the displacement is adjusted in this equation until tangency 
to the R-curve is achieved. In this way critical values are 
found for the displacement, S c , crack-extension, da c , and crack- 
extension force, G . The critical stress, a c , is found by 
Equation (8) or (3^, the load, P c , by Equation (5), and the 
toughness, K r , by Equation (1). The results for the first three 
tests are summarized in Table 2 . This table also includes 
measured values and results from the residual strength diagram 
discussed in the next section. There is good agreement. 


RESIDUAL STRENGTH DIAGRAM 


The residual strength diagram is used in an engineering approach 
to determine how the residual strength of a structure depends on 
the crack geometry and the specimen size. It is shown 
schematically for a center cracked sheet structure in Figure 6. 
The dotted line represents the residual strength assuming linear 
elastic fracture mechanics for an infinite sheet. The straight 
line shows the residual strength assuming net section collapse of 
a sheet of width W. For small and large crack sizes, near the 
intersections of these two curves there is a region of transition 
from one failure mode to the other. Feddersen [ 4 ] argued that 
two linear tangents to the idealized K curve can be used to 
establish a smooth and continuous curve for the residual 
strength. One tangent to the K curve is drawn from the point 
a=a where a y is the stress at which the structure without crack 
collapses. The other tangent is drawn from the point 2a-W that 
represents the width of the structure. The fracture toughness is 
customarily determined from the crack size at the failure stress 
by Equation (2) . However, in thin sheets where crack extension 
occurs before failure, the final crack size is generally not 
known. Therefore the data points plotted in the residual 
strength diagram usually are the initial crack size and the final 
failure stress, i.e. the crack extension and plasticity are 
ianored. To compensate for this inconsistency, an effective 
value of the toughness is used, denoted by K e , which is lower 
than the true toughness, such as that obtained e.g. by the R- 
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curve. A slight improvement to the analysis can also be made by 
taking the finite width of the panel into account, i.e. by using 

K e = a c t naL i sec(ira f /w) ]* ( 12 ) 

In this way a value of K e =lll MPaym was obtained as the average 
of the first three tests. Application of the residual strength 
diagram to the first three tests is shown in Figure 7. The 
predicted fracture ^ paths are shown and the measured data points 
are plotted as solid circles. The results for all three tests 
f a ll in the fracture mechanics region of the diagram. They are 
summarized in Table 2. There is good agreement between the 
measured data, the predictions from the R-curve and the residual 
strength diagram for these single crack specimens. 


THE LINK-UP CRITERION 


For the panels with multiple site damage (MSD) we have used the 
analysis of Broek [1]. Here we have a situation where small 
cracks exist ahead of the large central crack. The first link-up 
is governed by the stress-intensity factor of the combination of 
the main crack and the first MSD crack. The geometry factors of 
the two cracks . must therefore be compounded by the effect of 
crack interaction. As the ligament undoubtedly will fail by 
plastic collapse or net section yield, the compounded stress 
intensity factor can be used to calculate the plastic zone of 
both cracks and used with the criterion that link-up occurs when 
the two plastic zones meet. Broek found that the best estimate 
of the plastic zone size for this application was given by the 
Irwin expression 

r p = (K/ct y ) 2 /2tt (13) 

where <7 y is the collapse strength. This expression represents 
the plastic zone at the tip of the main crack and a similar 
expression holds for the plastic zone at the MSD crack. Here the 
stress-intensity factor is still given by Equation (2) , but with 
the MSD crack nearby the function 0 that takes the crack 
interaction into account is much more complicated than Equation 
(4) . To model the first link-up Broek assumed that the main 
interaction was between one tip of the central crack and the 
first MSD crack in front of it, ignoring all the other MSD 
cracks. . Thus will depend on the central crack size, a, the MSD 
crack size, a MSD , and the ligament L between them. Figure 2, so 
that Equation (2) becomes 

K = (3 (a, a MSD , L) oj(ira) ( 14 ) 
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A similar expression holds for the stress-intensity factor to use 
in the plastic zone at the MSD crack. For these fi functions 
Broek has derived approximations based on published expressions 
for the stress intensity factors of two unequal length collinear 
cracks, which we have used. The link-up criterion between the 
central and first MSD crack now is 


r p + r pMSD ^ 


or, from Equations (13) and (14) 

o — Oy [ 2 L/ ( a/3 + a MS0 ^HSo ) ^ 


(15) 


(16) 


For test MSD— 4 this expression is plotted as the first dashed 
curve in Figure 8. Though it is not shown explicitly in Equation 
(16) this curve represents the final failure stress as a function 
of initial crack size, in the same sense as was discussed for 
Equation (12) . The x-axis in Figure 8 represents the central 
crack size, a. For the curves that are plotted the MSD crack 
size, a MSD , and the position of the MSD crack, d MS0 , are kept 
constant. However, as may be clear from Figure 2, the ligament L 
varies with a. For the central crack with size a=177.8 mm. 
Equation (16) predicts the first link-up at a=81 MPa. This is 
shown by the intersection of the fracture path line with the 
first dashed curve in Figure 8. After the first link-up the 
central crack has absorbed the first MSD crack and it is assumed 
that we now have a larger central crack, a=195.58 mm, that 
interacts with the second MSD crack. Equation (16) is again used 
with these new parameters and plotted as the second dashed curve 
in Figure 8. The predicted stress for the second link-up then is 
a=117 MPa. This procedure is then repeated for the third link-up 
and we find a=lll MPa at a=220.98 mm. Since this stress is lower 
than for the second link-up, we conclude that the panel is super- 
critical and the fracture will run straight through after 
absorbing the second MSD crack at 117 MPa, as shown by the 
fracture path in Figure 8. When all three MSD cracks have been 
absorbed by the central crack, it behaves as a large single crack 
of size a=246 .38 mm in the panel. Failure is then predicted by 
the simpler Equation (12), as also shown in Figure 8. The 
results of this analysis for tests MSD-4, 5, and 7 are summarized 
in Table 3 and compared with the measured values. 


R-CURVE AND MSD 


We now apply the R-curve to the failure prediction of a panel 
with MSD cracks. We assume that an R-curve emanates from each 
MSD crack as well as the main central crack, as shown by the 
solid lines in Figure 9. Each R-curve is given by Equation (6) . 
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link h u»« d A«H ^ a " g . e 1 n .^ to each R-curve to determine the 

lmk-ups and final instability. The G-lines are still given by 

sl J llar to Equation (3), but now the p function for 
the central crack must also take the interaction with the MSD 
crack into account 


G(a) = (a 2 jra/E)/3 2 (a, a. 


MSD' 


L) 


(17) 


F °r a given stress this equation gives the G-line under load 
controi. Using the appropriate parameters, the stress was 
adjusted until tangency was achieved for each of the R-curves, as 
shown by the dashed lines in Figure 9. The final fracture at 

*5® 1 ® aine P rocedur e as discussed before. With the 
MSD cracks the G-lines have quite a bit of curvature just before 
the instability point. The link-up stress for the third MSD 
crack is less than for the second. This indicates that, after 
ii nk Z UP the .ff cond MSD crack, the panel is super-critical and 
M 9 n run straight through after absorbing the second 

in ; Sk res ? lts for tests MSD— 4 , 5, and 7 are summarized 

There good agreement with the measured values and 
also with the results predicted from the plastic zone model. 

the 6 n?«i-f^ r MSD " 5 2 n ? ? P red iction from both the R-curve and 
the plastic zone model is that the final fracture will occur with 

the second and third link-up, whereas actually the load had to be 

£f om 375 to 407 ^ for MSD— 5 and from 205 to 214 kN for 
MSD 7 after the second and third link-ups for the final fracture. 


CONCLUSION 


Analysis appears to provide good predictions of the residual 
strength, link-ups, and fracture of panels with multiple site 
damage (MSD) of different sifce and spacing. The critical 
fracture stress can be predicted by using either the plastic zone 

° r the R ” curve ' and the results are virtually the same, 
improvements to the analysis are being considered. 
Notwithstanding the large widths of the panels the MSD cracks are 
quite small and closely spaced so that small discrepancies in 
these dimensions may affect the results. Therefore more accurate 
measurements of the MSD configurations are currently being made 

a m i; ro f?°P e l T 5 e plastic zone analysis uses guestimates 
for the effective fracture toughness, K , and the collapse 
strength, a„. The values of these two quantities can be 
optimized when the data from all the tests are in. In the R- 

load control was used and the backward growth of 
the MSD crack was ignored. Improvements in this analysis are 
quite complex, but we think possible. Displacement control would 
allow the load to drop as in fact it did in MSD-5 and 7. 
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Table 1. Planned Test Matrix for 90 MSD Program 


Test 

No 

Main 

crack 

MSD cracks 



Date 

of 


2a 

(mm) 

a 

(mm) 

^MSD ®MSD 

(mm) (mm) 

2a MSD 

(mm) 

number 

per 

side 

test 


MSD-l 
MSD- 2 
MSD- 3 
MSD- 4 
MSD- 5 
MSD-6 a 
MSD- 7 
MSD— 8 b 
MSD- 9 
MSD- 10 
MSD-11 
MSD-12 


355.6 

177.8 



0 

203.2 

101.6 



0 

508.0 

254.0 



0 

355.6 

177.8 

190.5 

25.4 

10.16 3 

142.24 

71.12 

88.9 

38.1 

15.24 3 

508.0 

254.0 

(no anti-buckling guides) 

508 . 0 

254.0 

266.7 

38.1 

12.7 5 

355.6 

177.8 

190.5 

25.4 

10.16 3 

482 . 6 

241.3 

266.7 

38.1 

12.7 10 

254 . 0 

127.0 

165.1 

25.4 

10.16 10 


(to be defined) 
(to be defined) 


7/23/93 

8/25/93 

9/13/93 

11/29/93 

12/1/93 

4/12/94 

4/20/94 


* This test was recommended by NASA Langley 
This test is a repeat of MSD-4 


Table 2. 


Measured Data and Failure Predictions for the First 
Three 90 MSD Tests 


MSD-2 

Measured 

RSD 

R- curve 

MSD-l 

Measured 

RSD 

R-curve 

MSD- 3 

Measured 

RSD 

R-curve 


Half-Crack Size 


Stress Load Displacement 


Initial 

(mm) 

Final 

(mm) 

Extension 

(mm) 


^ HUll ) 

101.6 



184 

428 

16.1 


142 

40.4 

196 

454 



141.1 

39.5 

183 

424 

14.5 

177.8 



148 

343 

12.6 


244.3 

66.5 

146 

340 



247.1 

69.3 

146 

340 

11.9 

254 . 0 



124 

289 

12.2 


340.4 

86.4 

121 

280 



353.1 

99.1 

125 

289 

10.5 


Date 


8/25/93 


7/23/93 


9/13/93 
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Table 3 


Link-up and Failure Predictions from the Residual 
Strength Diagram and the R-curve for tests 90 MSD- 
4 , 5 , and 7 



Link- 

Load 

(kN) 

-up 1 
Stress 
(MPa) 

Link- 

Load 

(kN) 

-up 2 
Stress 
(MPa) 

Link- 

Load 

(kN) 

-up 3 
Stress 
(MPa) 

Fracture 
Load Stress 
(kN) (MPa) 

MSD-4 

Measured 

196 

84 

227 

98 

227 

98 

307 

132 

RSD 

188 

81 

272 

117 

(258) 

(111) 

285 

123 

R-curve 

240 

103 

274 

118 

(258) 

(111) 

294 

126 

MSD-5 

Measured 

320 

138 

375 

161 

375 

161 

407 

174 

RSD 

329 

142 

445 

191 

(392) 

(169) 

(346) 

(149) 

R-curve 

388 

167 

403 

174 

(347) 

(149) 

(344) 

(148) 

MSD-7 

Measured 

133 

57 

205 

88 

205 

88 

214 

92 

RSD 

134 

57.8 

268 

116 

(249) 

(107) 

(204) 

(88) 

R-curve 

183 

78.9 

255 

110 

(239) 

(103) 

(213) 

(92) 


The numbers in parentheses are lower than for a previous link up, 
thus failure will precede the link-up 
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Figure 1. Typical 90 MSD test panel 
2286 mm wide (90 in) 



Figure 2. Definitions of 90 MSD crack configurations and dimensions 
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G, R (kN/m) Load (kN) 
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Figure 3. Crack growth measurements for tests 90 MSD-1, 2, 3 and 6 
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Figure 4. Crack-extension force data and analytic R-curve for 90 MSD tests 





j Initial crack size: 2a 0 = 355.6 mm 

A 

Data Points 

▲ 

Fitted Points 


R-curvc: R = 75.9 da° 28 

— 

G c -linc. G(a) = *o c J a/E 


At Instability; 
a c = 146 MPa 
da c = 69.3 mm 
= 248 kN/m 


Crack Extension, da (mm) 
Figure 5. Stability diagram for test 90 MSD-1 


K e = 1 11 MPaVm 
cr Y = 345 MPa 
W = 2286 mm 


Net Section Collapse 


V(n a) 


Fracture 

Mechanics 


500 1000 1500 2000 W 2 

Crack Size, 2a (mm) 

Figure 6. Residual strength diagram with analysis of Feddersen 







Figure 7. Residual strength for tests MSD-1,2, and 3 with measured data points and predicted fracture paths 



Figure 8. Residual strength diagram for test MSD-4 with plastic zone criteria and predicted fracture path 
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G, R (kN/m) 








